Black Hole Production from High Energy Scattering in AdS/CFT 
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In this article we show how to set up initial states in TV = 4 SYM theory that correspond to high 
energy graviton collisions, leading to black hole formation in AdSs, x S . For this purpose, we study 
states in the gauge theory that are dual to graviton wavepackets localized at the center of AdSs, 
and carrying large angular momentum along the S 5 . These states are created by exciting only the 
s-wave mode of one of the complex adjoint scalars of SYM. For a single graviton, the state is 1/2 
BPS and one can show that it is dual to a linearized 1 /2 BPS geometry in the bulk. Exploiting this 
dictionary, we show how to localize the particle's wavefunciton so that the dual linearized metric has 
the form of a Aichelburg-Sexl shock wave. One can then put two such shock waves into a head-on 
collision, which is known to produce a trapped surface. Finally, we discuss the prospect of studying 
graviton scattering directly at strong coupling in the gauge theory using a reduced model of matrix 
quantum mechanics. 

PACS numbers: Valid PACS appear here 



INTRODUCTION 

Since the discovery of the AdS/CFT correspondence 
there have been great expectations that this dual- 
ity will help us answer deep questions about black hole 
physics. The reason is that the gauge theory provides a 
non-perturbative description of quantum gravity in the 
bulk with fixed boundary conditions. In particular, one 
expects to find a resolution to the black hole singularity, 
and a solution to the information loss problem. For the 
former, most research has been focused on finding signa- 
tures of the eternal black hole singularity in the gauge 
theory correlation functions [2| . This has been done ex- 
ploiting the duality between the gauge theory at finite 
temperature, and an eternal AdS-Schwarzchild black hole 
in the bulk. 

However, one would like to understand the dynami- 
cal phenomena of black hole formation, and its eventual 
evaporation. For this reason, it is very desirable to un- 
derstand how to set up initial states in the gauge theory 
which are dual black hole forming processes in the bulk. 
One way to do this is to study states in TV = 4 SYM 
theory which are dual to high energy graviton scattering 
in the bulk. 

Classically, and in flat space, it is known that for ener- 
gies much greater than the Planck's energy, these scatter- 
ing processes lead to the formation of a trapped surface 
0| . This has been shown by superposing two Aichelburg- 
Sexl shock waves, which represent the gravitational field 
of massless particles in flat space However, quan- 
tum mechanically one needs to be more careful and take 
into account the wave nature of the particle. This was 
discussed in Q. 

Scattering processes have also been studied in the con- 
text of the AdS/CFT duality in, for example 0]. How- 
ever, most of these works have dealt with the problem of 



particles coming from the boundary into the AdS bulk. 
This has various problems. The first is that the high 
redshift between the bulk and the boundary makes it 
hard to focus wavepackets on small scales. Secondly, 
since these calculation only involve correlation functions 
on the boundary, one is necessarily calculating S-matrix 
elements. It is desirable to have a precise state in SYM 
on R x S 3 , which we could evolve in time (at least in 
principle). Finally, focusing the scattering process into 
the center of AdS should lead to simplifications since the 
whole process is basically taking place in flat space. 

Therefore, in this paper we answer the following simple 
question: what is the SYM state dual to a Aichelburg- 
Sexl shock wave for a graviton well localized at the center 
of AdS, and traveling along the S 5 ? 

Moreover, we show how to superpose two such waves 
so we get a two-graviton scattering process in the bulk. 
We mainly work at zero impact parameter, but other 
generalizations are straightforward. We also make sure 
that the resulting trapped surface is much larger than 
the string scale, but much smaller that the AdS radius. 
Furthermore, we ensure that the curvature at the trapped 
surface is small. Finally, we comment on how one should 
understand this scattering process at strong coupling and 
at non-zero impact parameter, in terms of a quantum 
model of commuting matrices as in Q. 

The organization of the paper goes as follows. In sec- 
tion II, we take a 1/2 BPS linearized perturbation of 
AdS$ x S 5 described by the so-called LLM solution of 
9]. These solutions can be classified by "droplets" on a 
two dimensional plane. A droplet with the shape of a 
unit disk represents the AdS$ x S 5 solution. We then 
make a small ripple along the boundary of the unit disk. 
This represents a massless graviton wavefunction. We 
then take the limit of large radius of AdS while focus- 
ing on the region near the wavepackct. This leaves a ten 
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dimensional flat space metric with a linear perturbation 
of the Aichelburg-Sexl type. The difference here is that 
we can now control the spread of the wavepacket, and 
the usual delta function in the metric gets replaced with 
a smeared finite function. This allow us to control the 
curvature near the trapped which was an issue discussed 
in @. 

In section III, we write down states dual to these LLM 
shock waves using the dictionary recently developed in 
[irj| . The states are constructed out of the zero mode of 
one of the adjoint scalars of SYM theory. We give par- 
ticular examples of initial states corresponding to two 
gravitons separated by a finite distance in flat space, and 
heading towards each other with zero impact parameter. 
We also comment on the generalization for non-zero im- 
pact parameter. 

In section IV we speculate on how one should study 
these initial states and their time evolution, at strong 
coupling. In particular, we argue that the reduced matrix 
model of [8| should be sufficient to understand graviton 
scattering at impact parameters b > El 2 , where E is the 
center of mass energy of the collision. 

In section V we close with some final comments and 
future directions. 

Conventions. In this article we will set the ten dimen- 
sional Planck's constant l p — 1. In these units, the AdS 
radius is related to the rank of the SYM gauge group 
N by R = (471-iV) 1 / 4 . The string length is given by 
l s = \fa' = .g -1 / 4 , where g is the closed string coupling 
also related to the SYM coupling by A-xg — g\ M - We 
will also introduce the parameter h — X/N. 



The metrics can be written as, 



ds 2 = -h~ 2 (Dt) 2 + h 2 (dy 2 + dzdz) + yer G dtt 
+ye G dhl, 



= 2ycoshG, 
/ = -tanhG, 



f(z,z, y) 
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2 J 



d z z 



p(z',z') 
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Here, we have defined the covariant derivative, Dt = 
dt + V — dt+ ^iVdz — ^iVdz, and we are using complex 
coordinates in the y = plane: z = x\+ix%, z — x\ — ix2- 
Moreover, 



V(z,z,y) 



1 



dfz 



, p(z',z>)(z-z>) 
(\z - z'\ 2 + y 2 ) 2 



(2) 



All non-singular solutions must have p — ±l/ir. There- 
fore we can separate the integrations above in domains 
or "droplets" (T>i) for which p = 1/tt (say) inside and 
p = —l/n outside. In this paper we will only consider a 
single droplet at the center of the complex plane. 

The curvature radius of the azymptotic AdS?, x 5 5 (R) 
is set by i? 4 = J v d 2 z/n. Therefore, we can rescale all 
spatial coordinates by x l — » R 2 x l and we get an over- 
all factor of R 2 in front of the metric. Then, our area 
quantization condition is simply, 



(3) 



This makes it easier to compare with the gauge theory. 
Finally, the energy of the solution can be written as, 



THE AICHELBURG-SEXL METRIC FROM IIB 
SUPERGRAVITY 



E = J= tJ- / d 2 z\z\< 
n 2 TT J v 



1 



(4) 



In this section we will show how the Aichelburg-Sexl 
metric 0] can be derived from the 1/2 BPS geometries 
of type IIB supergravity @. This metric describes the 
external gravitational field of a point-like massless parti- 
cle. In fact, our version of the Aichelburg-Sexl metric will 
take into account the wave nature of the particle. This 
will allow us to control the curvature of the geometry 
near the region of the trapped surface. 

All 1/2 BPS solutions to type IIB supergravity with 
N units of Ramond-Ramond five-form flux have been 
found in They preserve a bosonic R x 50(4) x 50(4) 
symmetry of the ten dimensional space-time. Moreover, 
they preserve half of the supersymmetries of the original 
AdS§ x 5 5 . We will leave the details of this analysis to 
reference 0. Here we will only need the general form of 
the metric. 

It turns out that all solutions are classified by a single 
function, which we call p, on a two dimensional plane. 



Before turning our attention to the Aichelburg-Sexl 
metric, lets recover the familiar AdS 5 x 5 5 and flat space. 
For the former, we need to consider a circular droplet of 
unit radius. The integrations above give, 



y _ l + y 2 + \z\ 2 -^{l + y 2 + \z\ 2 ) 2 ~A\z\ 2 



Similarly, 



2V(l+2/ 2 +M 2 ) 2 -4|z| : 



2V(i + y 2 + M 2 ) 2 -4|z| s 



(6) 



Plugging these results back in the metric, and making 
the coordinate changes (z = re 11 ?), 



sinhpsin# , 
cosh p cos 9 , 
4>-t, 



(7) 
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we obtain the standard form for the AdSs x S 5 metric: 

2 2 2 2 ~ 2 

— cosh pdt + dp + sinh pdil 3 

~ 2 +am 2 ednf\ . (8) 



ds 2 = R 2 



+d9 2 + cos z 



Let us now obtain flat space. For this, we need to zoom 
in into the edge of the droplet. We will do it by rescaling 



(9) 

while taking the limit R — > oo. We can then make the 
following coordinate change: y — rir 2 , x 2 — (r\ — rf )/2. 
A careful expansion of the metric above gives, 



1 

R 



ds 2 



-dx\(2dt — dxi 



dr 



(10) 



where f = rj + r~2 is an eight dimensional transverse vec- 
tor to the droplet. Note the difference between this limit 
and the usual pp-wave limit. In the later, one does not 
rescale the time coordinate and, instead, make a further 
rescaling x\ — > x%/R. 

From the metric (fT0| we clearly see the two null 
geodesies at the center of AdS: x\ = const., and x\ 
'21 + const. Ripples on the boundary of the droplet fol- 
lowing any of these geodesies will be 1/2 BPS, but with 
opposite angular momentum on the S 5 . Let us now focus 
on a ripple centered at Xi — 0. We can later boost this 
result to obtain a particle along the other geodesic. 

The ripple will be given by a fluctuation on the radius 
of the droplet: 

^boundary = 1 + Sr((f>) . (11) 

The area quantization condition |3]) implies, 



d0<5r(</>) + / dcj> [6r((f>)] 



. 



(12) 



To leading order in the perturbation, one can set the 
second term above to zero in most calculations. However, 
for the energy ((4]) one needs to be careful. Using (p~2|) in 
Q we get, 



(13) 



Now we would like to localize the perturbation Sr 
within the flat space patch near <fi — tt/2 (or X\ = 0). 
Moreover, one would like to keep the energy E fixed as 
we take the R — -> oo limit. To this end, we will normalize 
the solution as 



5r((f>) = 



j^g[RU- 

=— g(xi/X) . 

i? 3 VA 



-/2)/A] 



(14) 



Here the function g (x) is zero outside the interval x S 
[—1/2, 1/2], and A is some wavelength that does not scale 
with R. Moreover, the function g(x) is normalized such 
that, 



1/2 



1/2 



dx[g(x)} 2 = 1 



Using the area constraint (J2|) , we also get 



1/2 



dxg(x) 



-1/2 



2tt 3 / 2 V£ 
i? 3 VA 



0, 



(15) 



(16) 



where in the last step we have taken the large R limit. 

It is easy to verify from (fl"3|) that this perturbation 
indeed has energy E, once we take into account the fact 
that we have rescaled the time as in ([9]), and hence the 
Hamiltonian as H ncw = H \d/ R- 

Let us now derive the metric corresponding to this per- 
turbation, starting with the function / of the metric ©. 
We can write / = /o + jj?5f, where /o is the AdS space 
result Eq. ©. Similarly, V = V + SV, where V is 
also the AdS space result ©. The reason for the scal- 
ings with R will become apparent in a moment. For now, 
we just expand the metric to linear order in the pertur- 
bations and take the R — > oo limit with the change of 
coordinates (GO). The result is, 



ds 2 = ~dxi(2dt — dx\) + dr 



7?2 



Sgudxf , (17) 



where 



Sgn = i(6V - SV) 



2x2 Af 

y 



(18) 



Let us now study the form of these perturbations more 
closely. One can easily show that, after doing the rescal- 
ings ([9]), the metric perturbation becomes: 



Sg n = 



2R 3 


r \/2 


fR 2 Sr(x' 1 




-J 


dx\ 






7T 




-A/2 


L _ 




X 




+ |«- 


a?) 




x 


y 


f R 2 (x 1 - 


-x\) 


2 + 



dXn 



{X2 - X ' 2 -^ X 1- X f)f 



0(1/ R)) 



(19) 



We can simplify this expression in two ways. First let 
us note the identity 



R 



R™oo (1 + [RG(x)} 2 



R 



de i?™<L e + [RG{x)] 2 
d_ J_ 
"de 7- 



,ttS(G(x)) 



-S(G(x)) 



e=l 



(20) 
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Moreover, from (fl4|) we can see that R 2 Sr ~ 1/R and 
so it is very small for x 2 ~ 0(1). Therefore, using ([M]) 
and making the usual change of coordinates y — rir 2 , 
x 2 = 0? - r|)/2, we get 



$9 



11 



(4tt) 3 £; 



5\(xi) , 



where we have defined, 

,1/2 



8\ (x) = I dx l [g(x')] 2 S(x — \x) 
-1/2 



(21) 



(22) 



This is basically a delta function with support in an in- 
terval Ax = A. 

To summarize, we have found that for a perturbation 
of the unit disk of the form (fT4|) . the LLM metric in the 
large R limit becomes, 



-dx\{2dt — dx\ 



dr^+ i -^^6,(x 1 )dxl 

n 



(23) 

This is precisely the form of the Aichelburg-Sexl metric in 
ten dimensions :3j . Note that the perturbation obeys the 
ten dimensional Laplace equation with a delta function 
source at r = 0. This delta function is only an artifact 
of the limit R — > 00. The full ten dimensional solution is 
completely smooth even at r = 0. 

We can then see that for a 1/2 BPS perturbation with 
the opposite angular momentum one can simply shift this 
result by x\ — ► It — x\. At the level of linearized grav- 
ity, we can superimpose both waves and get a solution 
to Einstein's equations. Of course, the solution will not 
be BPS anymore. Moreover, the simple superposition is 
only valid outside the region X\, 2t — x\ > 0, where the 
colliding shocks start to influence each other. The full 
solution inside this region is still not known. 

The two shock waves will collide at t — and at x\ = 0. 
It is known that such head-on collisions lead to the for- 
mation of an eight dimensional closed marginally trapped 
surface with radius [3j 



(24) 



We now need to ensure that this classical calculation is 
not spoiled by high curvatures near the trapped surface. 
This issue was discussed in jfj. For a single shock wave, 
the only non-zero components of the curvature tensor are 



I d 2 (AnfE 

R ^ = -2 5x{Xl) dnd rj |rf • (25) 



All curvature invariants are finite unless at the precise 
location of the particle r — 0. 

For the two-wave superposition, one encounters a sin- 
gularity at t = 0, xx = if we take the limit A — > 0. In 



fact, in this limit the only diverging curvature invariant 
is of the form 191 



(26) 



where we have taken into account the peak in the function 
8\(x). We can evaluate this quantity at the radius of 
the trapped surface (|24|) . and require: (i? AiypCT ) 2 <C 
\fh\ ^> A and \ fh\ 3> l s - These conditions translate into 

. g-WE-W « A « E 1 ' 7 , 

• E 1 ' 7 » g' 1 ^. 

It is clear that for large energies these constraints are 
very easy to satisfy and they allow us to take A <§; 1. 

At this point, we would like to emphasize that in order 
to obtain the flat space limit with a fixed string length, we 
cannot take the usual 't-Hooft limit. That is, we do not 
take gy M N = fixed. Instead, we take N — > 00 keeping g s 
small but fixed. Indeed, this is basically the same limit 
one takes when studying pp- waves [l2| . 

INITIAL STATES FROM Af = 4 SYM 

In this section we will review the basic dictionary be- 
tween the 1/2 BPS states of TV = 4 SYM theory on 
IxS 3 , and matrix eigenvalue distributions. These eigen- 
value distributions can then be mapped to the LLM ge- 
ometries studied in the previous section [9(. Using this 
map it is easy to find the dual state to the Aichelburg- 
Sexl metric (|23|) . Moreover, it is straightforward to write 
a superposition of two shock waves. 

The 1/2 BPS states of Af = 4 SYM theory are created 
with the zero mode of one of the complex adjoin scalar 
field Z(t, O), where fl £ S 3 . The tree- level quadratic 
action for this mode is simply 



S = / dt Tr 



Z\ 2 - 



(27) 



This simple model has been discussed in [13|. 

We can now introduce the usual oscillators. We start 
with the following canonical commutation relations, 

\{Z)l (IT)'] = i5\5{ , [{Z)i (U) l k ] = l S l l S{ , (28) 

where EE = Z, and II = Z. The harmonic oscillator 
operators are defined as usual, 



A»--jj(*-ffl) 



i* = i(2-in) 



with 

Note that IT = LT 



k] = 8 l l 5{, [Ai,(£) 



l k ] = s\si, 



(29) 



(30) 
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The Hamiltonian of the system is simply [26J , 

H= Tr(AU + IU) . (31) 

Note that these operators carry a single unit of the 
U(l) C SU(4) R-charge. The generator for this U(l) 
is 



J = Tr(A*A-A*A) 



(32) 



So we see that A^ has charge +1 and A^ has charge —1. 
This R-charge is dual to an angular momentum along the 
S 5 in the bulk. Since [H, J] — 0, one can define a new 
Hamiltonian, 



H' = H — J 



(33) 



which is naturally identified with the time t in the LLM 
coordinates used in the previous section. 

The (anti) 1/2 BPS states are formed by acting with 
only (A^) A^ on the ground state. In the time slicing ([33]) . 
the 1/2 BPS states are time independent and the anti-1/2 
BPS have a simple time dependence where A^ — > e 2lt A^. 

Lets start by studying "coherent states" of the form, 



1^1/ 



2 BPS 



) oc e 



Trn(At) 



|0) 



(34) 



We can also build non-BPS states that include two wave- 
functions with opposite R-charge: 

l^oce^^e^^lO) . (35) 

In the dual string theory, these states will correspond 
to two gravitons at the center of AdS$ and traveling in 
opposite directions along the equator of S 5 . 

Although the state ([33)1 is not BPS, as long as the 
gravitons are initially well separated in the S 5 , one can 
treat each exponential as a separate initial 1/2 BPS state. 
This will allow us to use the dictionary between the 1/2 
BPS states (|34| and the LLM geometries developed in 
(Io| . Let us now review this dictionary. 

We can start by defining coherent states for the an- 
nihilation operator A. Namely, A\\Z) = Z\\Z) and 
where the resolution of the identity takes the form 1 = 
J[d 2 Z]\Z)(Z\, with [d 2 Z] the U(N) invariant measure 
over complex matrices. 

The 1/2 BPS state ([34]) becomes, 



Trfi(z)/ft e - Tr\z\ 2 /(2h) = MZ)e~ Tl \ z \ 2 /( 2h ) , 

(36) 

where we have rescaled the Z matrix for later conve- 
nience. By going to an eigenvalue basis, one can show 
that the normalization of this state is given by the par- 
tition function [lfj, 



N 



(37) 



where 



W(z, 



(38) 



It is well known that in the "classical" limit H — > 0, 
one can use the saddle point approximation and replace 
the sums in (|37[) by integrals over an eigenvalue density 
in the complex plane [l5| . This is the usual 2D Coulomb 
gas problem. 

The "free energy" functional is, 



F[p] = -~ J d 2 zp(z)W(z,z) 



d 2 z / d 2 z'p{z)p{z')\og\z- z'\ 2 .(39) 



1 

~2 

The saddle point equations are SF[p] — subject to the 
constraint J d 2 zp — fir 1 . Using the fact that log|z| 2 is 
the Green's function in two dimensions, one finds that 
the density p is locally constant with p = I /(fin). The 
eigenvalues are then distributed over these constant den- 
sity droplets which can be translated to the LLM geome- 
tries. Compelling evidence for this was specially given in 
IcL 3] , where it was shown that one can recover a probe 
string sigma model on a LLM geometry in the gauge the- 
ory by including a single trace "probe" state in top of the 
1/2 BPS state flM|)- 

The energy of a 1/2 BPS state can be calculated by 
the free field theory Hamiltonian (pH]) . Using the coher- 
ent states, the Hamiltonian is simply H — Tr(ZdZ), 
and it acts on the wavefunction ip(Z). Doing a partial 
integration in the Random Matrix Model, it is easy to 
show that the energy of the state is given by, 



E 



Tr\Z\ 



1 

2fi2 



1 

fi 2 ^ 



d 2 z\z\ 2 



v 



1 



(40) 



This expression matches precisely with the energy calcu- 
lated in the gravity side, Eq. 

We can now take a potential of the form 



n(z) 



k>0 



(41) 



It is known that for a single droplet the moments are 
related to the shape of the droplet by [HI 



1 



dz 



(42) 



For a linearized perturbation in the boundary of the form 
(fTTjl . the moments reduce to, 



t fc » -L d4>8r(4>)e~ lk * 

7Tfc ,/n 



(43) 



Therefore, we can find the gauge theory dual to the 
Aichelburg-Sexl metric lp3]) . We simply use (H} in (|43jl : 

tk « L[B: e -**/* f 1/2 dxg{xy kxX ' R . 
k V 7T J_ 1/2 



(44) 
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The single particle state can now be written as, 



exp 



fl/2 , 

/ dxg(x) Tr log ( ie~ ixX > R - VhA^ 

J-l/2 V 



|0). 

(45) 



It is now very easy to write down the anti-BPS state 
which is located at the null geodesic x\ = L at t = 0. We 
simply shift x — > x + L and replace A** — > A^ inside the 
log of (|45p. The two particle initial state is simply the 
product of two such functions acting on the ground state 
as in (l35l). 



The Problem with Single Trace States 

In this section we would like to point out an obstruc- 
tion in constructing the Aichelburg-Sexl metric using 
only single trace states in SYM. 

Consider the 1/2 BPS state 



n>0 V 



(46) 



In the large N limit, this state has unit norm so long as, 

E i^i 2 = 1 ■ ( 4? ) 

n>0 

Of course, this is only valid as long as the amplitudes 
\4>n\ 2 decay sufficiently fast for n > n max _, where n max <C 
N. This constraint comes from the fact that, for suffi- 
ciently large n, the single grace states with different ns 
are no longer orthogonal. 

The "free energy" functional that comes from comput- 
ing the square amplitude using coherent states is 



i 



d 2 zp(z)\z\ 



-~ / d 2 z I d 2 z'p(z)p(z')log\z-z'\ 2 



log 



J \n>0 v 



(48) 



The saddle point equation <5-F[p] = with the constraint 
J d 2 zp — ftr 1 gives, 



const. = — — —\z\ 
h \ 



J2tkZ k + h.c. j 

fc>o / 



d 2 z'p(z')log\z- 



(49) 



where the constant comes from the Lagrange multiplier 
implementing the area constraint, and the moments tk 
are given by, 



tk — 



Vk(fd 2 zp{z)J2 n> o^=< 



Vfc (E n>0 Vntpnt- 



(50) 



In the last line, we have used the definition of the mo- 
ments in terms of the eigenvalue distribution, Eq. (|42jl . 
Moreover, at the linearized level, one can see from (|43|) 
that t- n « -t*. 

One can then use ([50]) and the normalization (|4T|) to 
solve for the wavefunction: 



i\)k = 1 Vkt k , 

n 2 = E fc i^i 2 - 



(51) 
(52) 



n>0 



The last equality comes from the normalization of the 
wavefunction, and gives an extra constraint on the mo- 
ments of the distribution. Note that we did not have this 
constraint with the exponential states. This new con- 
straint is part of the problem. To see this, lets construct 
a particular example. 

Consider a perturbation 8r(<f>) of the form (fl4|) with 



1 x G [—1/2, 0) 
j(x) oc { -1 x e [0,1/2] 
otherwise 



(53) 



Using the linearized form of the moments (|44j) one can 
easily finds, 



AihR E 



k 3 / 2 V ttA 



e '^/ 2 sin 2 



ARJ 



(54) 



The normalization of the wavefunction gives a relation 
between E and A: 



E 



A log 2 



(55) 



To obtain this relation, one approximates the sum over 
k in ([52"]) by an integral. 

The relation (|55[) is not present for the exponential 
states studied above. Moreover, we see that this is not 
compatible with the constraints studied in the previous 
section. In particular, we want to take E 3> 1 to have a 
macroscopic trapped surface. This means that we need 
A <C 1. In this case, however, the low curvature condi- 
tions would imply A 6 / 7 3> g^ 1 ^ 2 which is not compat- 
ible with l s /lp ~ 5~ 1//4 ^> 1- Therefore, we conclude 
that single trace states are not good enough to produce 
macroscopic black holes. 

There is another problem with single trace states: the 
ten dimensional wavefunction on the bulk is not well lo- 
calized in the orthogonal directions to the shock wave. 
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We show this in detail in the Appendix, by relating the 
single trace excitations to single particle states in the 
bulk. However, one would like to understand this prob- 
lem directly from the gauge theory. So far we haven't 
found a clear answer to this question. What seems to 
be happening is that the saddle point approximation is 
not good when we use the single trace states. That is, 
due to the logarimic term in the "free energy" (|48|l . the 
eigenvalue distribution has a lot of variance around the 
saddle point. It would be interesting to show this by a 
numerical simulation as in [16j |. 

Nevertheless, we know that this does not happen for 
the exponential states ([34]) . Indeed, we know how to com- 
pute the corrections to the saddle point in this case. This 
is discussed in [15| . The leading correction is a smooth- 



ing of the density distribution p(z) near the boundary of 
the droplet. That is, the density is not really a step func- 
tion, but it has an exponential decay at the edge of the 
droplet over a distance of order ~ \f% (in units where the 
circular droplet has unit radius). This correction is even 
present for the ground state circular droplet. If we trans- 
late this to the flat space coordinates (O, this means that 
the variation of the density occurs on scales of a; 2 ~ 0(1). 

Why don't we see this effect in String Theory? The 
reason is that, the string scale is very big l s = g^ 1 ^ 4 
compared to l p = 1, and so the string sees a step function 
to a very good accuracy. Therefore, the supergravity 
solution captures only the average shape of the droplet. 
This applies to general droplets. In other words, we only 
care about the shape of the droplet and not its small scale 
structure along the edge. The shape of the droplet will be 
accurately captured by the saddle point approximation as 
long as the variation of the shape of the edge occurs on 
a scale Ac/) ^> \fh. That is, we need a "long" wavelength 
variation. In our case, we saw in the previous section that 
A<p ~ 1/R ~ h 1 ' 4 > Vh. Therefore, we are allowed to 
use the saddle point approximation for the exponential 
states. 



THE STRONG COUPLING DESCRIPTION 

In the previous section, we have seen how to construct 
initial states in SYM theory which are dual to high en- 
ergy graviton collisions on the S 5 . However, this has been 
done using only the weak field theory description. For the 
case of states with the same R-charge, this description is 
enough since they are supersymmetric. However, when it 
comes scattering processes, one is dealing with non-BPS 
states and therefore one needs to work directly at strong 
coupling. In this section we argue that graviton scatter- 
ing at non-zero impact parameter should be described 
by a reduced quantum mechanical model of commuting 
matrices along the lines of 0]. 

Lets begin by reviewing the proposal of [8j . One starts 
by studying the dynamics of the zero modes of the six 



hermitean scalars of SYM, X a (a = 1,2,..., 6). The 
classical hamiltonian for these modes is given by [l6j |. 

H eff = T?(j2\{D t X a Y + l -{X a Y 

\a=l 



^f[X a ,X b }[X a 7 X b ] I (56) 



,,6=1 



8tt 



On general grounds, one expects the size of the SYM 
ground state to be, 



1 6 

-]T(0|Tr(A a A a )|0)~A 



(57) 



a=l 



in the large N limit. Therefore, if we re-scale the scalar 
by X a — > ^/NX a , we see that the commutator term in 
(|56|) will have an extra gyu^ factor in front of it. It then 
makes sense that, in the large N limit (or in the large 't- 
Hooft limit), the lowest energy fluctuations around the 
ground state should be described by a quantum mechan- 
ical system of commuting matrices. 

For commuting matrices we can go to an eigenvalue 
basis. In this basis we can define a vector containing the 
eigenvalues of all the matrices, afj = (A*-, . . . , A^), where 
i = 1, ...,N. The quantum Hamiltonian that follows 
from (1561 in this basis becomes [811 , 



i\ 2 1 ■ 



where 



M 2 =n 



Xi Xj I , 



(58) 



(59) 



i<3 



comes from the Jacobian produced from going to the 
eigenvalue basis. 

The wavefunctions must be symmetric under the ex- 
change of the eigenvalues by gauge invariance. Moreover, 
the measure in this basis becomes, 



(60) 



We can get rid of the extra factor of fi 2 in the inner 
product by rescaling the wavefunction as ip — ► ip/ fJ>. After 
this rescaling, the Hamiltonian becomes, 

ff = E(-^v,i + i|^) 
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where, 



6 E \Xi -X.j\ 2 

\%i ~ x j) ' { x i ~ x k) 



i j,k^i 



X% x j\^\ x i x k\^ 



(62) 



The ground state of this Hamiltonian is exactly known 
and takes the form [8(1 



■00 ~ /i ex P 



■5EI 



(63) 



In the large N limit, this wavefunction leads to an eigen- 



value distribution which consist of an S C 



The 



radius of the sphere was calculated in [l7| and it is given 
by = \J N/2. This confirms the expectation (|57|) . 

Now we want to add the dynamics of the off-diagonal 
elements of the matrices, for both the zero modes, and 
the higher spherical harmonics on S 3 . For this we fol- 
low the prescription of 0, [13, 0, 0, H3| , and quantize 
the off-diagonal modes in the background of the diagonal 
eigenvalue distribution. More specifically, one can derive 
an effective quadratic Hamiltonian for these modes that 
looks like 



HoS d 



lag. 



EE 



wtHAi)l{A a )) 



1,3 



with 



(«) 
y- • 



9ym I 
2vr 2 1 



(64) 



(65) 



The operators A a , A' a obey the usual commutator rela- 
tions for harmonic oscillators. The collective index a 
goes over the different kinds of excitations (scalars of 
fermions), and their spherical harmonic numbers. Fur- 
thermore, m a is their tree level mass, which is a posi- 
tive (half) integer for (fermions) scalars. The coupling to 
the eigenvalues comes from the commutator terms in the 
SYM lagrangian. 

Of course, with no supersymmetry one expects that 
integrating out the off-diagonal modes will produce a 
non-trivial functional determinant that will badly modify 
the Hamiltonian of the eigenvalues. The assumption of 
[H is that, given all the supersymmetry of SYM theory, 
the functional determinants between bosons and fermions 
should roughly cancel to leading order. It would be in- 
teresting to prove this rigorously. 

" dm 



Now, the proposal of 17 



is that we should 
treat the random numbers given by the distribution 
of the ground state (|63|) . Given the size of he ground 
state tq = ■J N/2, one can see from ([54]) that the off- 
diagonal modes will generically have masses of the order 



"T-off-diag. ~ \J 9ymN, and will decouple in the large 't- 
Hooft limit, or large N limit. 

It was first found in [l7| that, using this approach, one 
could reproduce the energy of the so-called BMN states 
to all loops in the 't-Hooft coupling. Moreover, a picture 
of semi-clasical "string bits" emerged which was later on 
confirmed in the classical string theory by the so-called 
"Giant Magnons" of Maldacena and Hofman 2l| . 

Now, let us suppose that we perturb the ground state 
by creating two waves of the form, 



(66) 



Let us also suppose that each individual wave fi a is "al- 
most" BPS. By this we mean that, if we had only one of 
them, the state would be very close to a BPS state. For 
example, one can show [l6| that holomorphic functions 
like fl ~ J2k tk zk are approximately BPS in the sense 
that (H — J) w in the large N limit. 

Finally, lets suppose that the two waves create two 
density fluctuations in the eigenvalues 5p a {x) that are 
both localized within a "flat space" patch. If we normal- 
ize the radius of the S 5 eigenvalue distribution to one, we 
have seen in section II that the flat space patch is located 
within an angular difference of A<j) ~ 1/R ~ iV -1 / 4 . The 
two density fluctuations will be located within this patch. 
That is, their angular separation is ~ b/N 1 / 4 , where b is 
a fixed number. 

We can now ask: are there any light off-diagonal modes 
that can be excited? One can divide the off-diagonal ex- 
citations into two categories. First, there are the "string 
bits" that connect different eigenvalues within the same 
density fluctuation Sp a . These should be suppressed since 
we assumed that each individual wave is almost BPS. 
However, we also have the string bits that join eigenval- 
ues between the two different lumps. From the disper- 
sion relation (l65|) , it is easy to see that these off-diagonal 
modes will have energies of the order 



E, 



off-di; 



ATi/4 



(67) 



However, we need to remember, that we are looking at 
small time scales and that we re-defined our time coordi- 
nate by t — > t/R, and so the Hamiltonian by H — > HR. 
Therefore, the energy of the off-diagonal modes in the 
new time coordinate is i? ff-diag > \fg~b = lj 2 b. To be 
able to ignore these modes, one would need to have large 
impact parameters so that E < -E ff-diag, where E is the 
center of mass energy of the collision. This leads to the 
bound, 



b>l 2 a E 



(68) 



In the string theory picture, this is precisely the bound 
to create long strings [22j . It is interesting to reproduce 
this bound from a gauge theory perspective. 
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However, to create black holes one needs b ~ ~ 
i? 1 / 7 , which is much smaller than (|68|) . This seems to 
suggest that the matrix model can only capture graviton 
scattering in the Born approximation. Nevetheless, there 
is evidence that stringy physics does not play a role in the 
formation of the black hole 22, 23[. It might be that the 



matrix model is good even for impact parameters lower 
than ([68]) . In any case, it would be interesting to even 
recover the 1/r 2 gravitational potential from the matrix 
model. 

Now let us explain in a bit more detail, how one should 
attack these problems. First, of all we need to learn how 
to write down well localized initial states using the SO(6) 
matrix model. The guess would be that one can take the 
same fi that we used for the single matrix model in sec- 
tion III. This needs to be checked in more detail. In 
particular, one would like to see if the resulting density 
fluctuation (for a single graviton) is indeed localized as 
desired. This might involve the need of numerical simu- 
lations as in 

To control the impact parameter using our formalism, 
one can start with an holomorphic wave of the form 
il ~ tkZ k (and the corresponding anti-holomorphic 
one), and perform a rotation of the coordinate z. This 
will take one of the geodesies off the equator. Then, 
by controlling the initial angle from which we start our 
lumps, one can arrange an initial state with two lumps 
that, after time evolving, will miss each other by a finite 
impact parameter. 

Finally, one needs to solve the time evolution of the 
density lumps Sp a using the Hamiltonian Q61[) . Note that 
this is a non-trivial A-body problem. However, one ex- 
pects some simplifications in taking the large N limit. 



CONCLUSION 

In this paper we have made a first attempt to care- 
fully define initial states in SYM theory that, under time 
evolution, are dual to high energy graviton scattering in 
the bulk. The basic approach has been to concentrate 
on the zero modes of the SYM scalars, which are dual to 
excitations that are located at the center of AdSs, but 
carrying angular momentum along the S 5 . We have given 
a dictionary between gauge theory states and linearized 
LLM geometries. Furthermore, we have shown how to 
obtain the Aichelburg-Sexl metric in the bulk, and how 
to superpose two such shock wave to produce a classical 
trapped surface. 

We have also paid special attention to suppressing 
stringy effects in the bulk, by producing a smoothed- 
out version of the shock- wave geometry. Finally, we have 
given evidence that scattering of gravitons a finite im- 
pact parameter is described by a reduced model of matrix 
quantum mechanics. 

There are many interesting question one can try to 



answer along these lines. First, one would like to learn 
how to construct the bulk wavefunction dual to the ex- 
ponential BPS states For single trace states, the 
dictionary is discussed in the Appendix. However, for 
these states, we showed that the wavefunction is not well 
localized in the directions orthogonal to the propagation 
of the particle. We argued that this is not a problem for 
the exponential states based on the validity of the sad- 
dle point approximation in the Random Matrix model. 
Nevertheless, it is very desirable to have an actual wave- 
function in the bulk to make an explicit calculation. 

Another interesting avenue of research is the study 
of gravitons scattering using the reduced matrix model 
studied in the previous section. In particular one would 
like to understand how to time evolve density fluctua- 
tions dp in the eigenvalue distribution. According to the 
arguments given above, we expect that such scattering 
processes will capture the physics of graviton scattering 
at large impact parameter in flat space. It would be very 
interesting to confirm this. 

The reduced matrix model could also help to answer 
the question about the wavefunction localization. For ex- 
ample, one can study the a holomorphic wavepacket of 
the form exp(^ ■ J2k and use the full SO(6) eigen- 

value measure to study its spread along all the directions 
on the S 5 . This will require the use of numerical tech- 
niques such as in 

It would be very interesting to understand how to cal- 
culate the finite time propagator (-01 (i) |-02 (0)} in the 
matrix model. In fact, we expect that for short times 
t ~ and for density fluctuations localized within an 
angular distance A(f> ~ this computation will give 

us the flat space S-matrix. One can then start to look 
for signatures of the black hole formation. In particular, 
one expects that these amplitudes should behave like 



(Mt)\M0)) 



(69) 



for large energies E [24[. Finally, note that the small 
black hole has a finite entropy Sbh ~ E s ^ 7 in the large 
N limit. Therefore, the matrix model (even without the 
off-diagonal excitations) has more than enough degrees 
of freedom to account for its microstates (I thank Rob 
Myers for pointing this out). 
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APPENDIX: THE TEN DIMENSIONAL 
WAVEFUNCTION 

In this appendix we study the relation between the sin- 
gle trace wavefunction (|4"6")l in SYM, and the ten dimen- 
sional wavefunction in the bulk. The main conclusion 
is that the single trace states lead to wavepackets that 
are not well localized in the direction orthogonal to the 
propagation of the graviton. 

To find the ten dimensional wavefunction, we need to 
linearize the LLM solutions, and quantize the fluctuation 
around AdS§ x S 5 . Fortunately, this was done in [25[. 
Their paramctrization of the edge fluctuation Sr((f>) takes 
the form, 



so that 



H4>) = J2 



(70) 



where a_ n = a* . 



After quantization, one finds that the Fourier coeffi- 
cients a n obey the commutation relations (after proper 
normalization) , 



[a n ,a) m ] = nS ntm . (71) 
This is very similar to the expectation value, 



(0|ft™ /2 r TrA n h m ^ Tr^nO) = nS n>m . (72) 
Therefore we can identify, 

4|0)~ Tr(At)»|0). (73) 
The metric fluctuation in global coordinates looks like 



S 9 ~ X c "^" s r 



where c n are constants, and 
Y n = e m ^-*) cosl"! 6 , a. 



cosh 1 " 1 p 



-a n e 



(74) 



(75) 



The important point here is that these modes satisfy the 
ten dimensional Laplace equation: 



v^ (s„r„) = 0. 



(76) 



Therefore, they provide a position representation of the 
wavefunction corresponding to the mode aJjO). In other 
words, 



(zioKJO) ~ s n Y n 



(77) 



Using the identification (|73|) , we see that the modes 
s n , Y n also give ten dimensional interpretation of the sin- 
gle trace states in the gauge theory. 

Let us now define the normalized modes, 



(j) n {x) (X S n Y n 



(78) 



(79) 



Here, E is a constant time surface. More explicitly, 



v /(n + l)(n + 2)(2n-l)(2n-3)cos"g in$ 
0n cosh"p 

. (80) 

There is a subtlety with the n = 1 mode which is not 
normalizable. In the LLM picture, this mode corresponds 
to a translation of the droplet in in the complex plane. 
We can set this mode to zero by translating the droplet. 
For example, shifting the whole droplet by a complex 
number a, the moments change by 



t fe -> t fe + a* 4,1 - a(k + l)t k+1 + 0{a 2 ) 



(81) 



Since the moments that we are considering are very small, 
we can just set a* — —t\ and get rid of this mode. The 
other modes will not be changed to leading order. 
A general normalized wavefunction is, 



n>l 



with 



n>l 



(82) 



(83) 



Given the relation between the bulk and gauge theory 
states (|73jl. it is natural to identify ip n above with the 
corresponding coefficients used in the single trace gauge 
theory state (|4"6"]) . Let us now use the example studied in 
the previous section. 

To study the spread of the wavefunction, is is conve- 
nient to make some coordinate transformations. Using 
the coordinate change (JT]), the rescalings © and the re- 
lation cos = r one can easily show, 



cos 8 
coshp 



1 - 



2R 2 ' 



(84) 



where r is the eight-dimensional vector transverse to the 
propagation of the graviton. Therefore, the probability 
density in a spatial slice S is, 



E V* V( n + + 2 )( 2 ™ " l)(2n - 3) 



fe>i 



x 1 



' e in<t> 



2R 2 J 



(85) 



Here, we have included the measure for the eight trans- 
verse directions. 

We can now look at the example studied in section 
III. Namely, the wavefunction given by Eq. (j5"4")l . After 
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integrating out the <j> coordinate, one finds that the full 
wavefunction (|85|) is peaked around \r\ ~ yR, with a 
spread of similar size. This is way too big! Therefore, we 
conclude that this kind of wavefunction cannot localize 
the graviton within the flat space region. 
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